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> ! ABSTRACT 
O 

^ ; We perform an adjust to the most recent structure function data, considering the QCD dipole 
picture applied to the ep scattering. The structure function F 2 at small x and intermediate 
Q 2 can be described by the model containing an economical number of free-parameters, which 
encodes the hard Pomeron physics. The structure function Fl and the gluon distribution 
are predicted without further adjustments. The data description is remarkably good, however 



■^j- ' giving an either low value for the fixed strong coupling constant. This indicates that a resummed 
next-to-leading level analysis should be done, which would bring it to suitable values. 

O 

Oh: 

j§ ■ 1 Introduction 



Deep inelastic electron-proton scattering experiments at HERA have provided measurements 



of the inclusive structure function F 2 (x, Q 2 ) in very small values of the Bjorken variable x (10 2 



down to 10 -5 ). In these processes the proton target is analyzed by a hard probe with virtuality 
Q 2 = —q 2 , where x ~ Q 2 /2p.q and p, q are the four-momenta of the incoming proton and 
the virtual photon probe. In that kinematical region, the gluon is the leading parton driving 
the small x behavior of the deep inelastic observables. The small x region is described in a 
formal way using the summation of gluon ladder diagrams, whose virtual contributions lead to 
the gluon reggeization. Such ladders diagrams are associated with the Pomeron, the leading 
Reggeon with the vacuum quantum numbers, which was introduced phenomenologically to 
describe the high energy behavior of the total and elastic cross-sections of the hadronic reactions 
and connected with the existence of large rapidity gaps in the produced final state [1]. 

At the leading order (LO) all powers of a s \ia(Q 2 / fj, 2 ), with [i 2 the factorization scale, are 
summed by the DGLAP evolution equations [2], which take into account only the strongly 
ordered parton transverse momenta (kx) ladders. At the present, the next-to-leading order 
(NLO) contribution is also considered, including non-ordered kx contributions in a covariant 
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way. In the current HERA kinematical regime the DGLAP approach is quite successful, al- 
though the theoretical expectation of deviations due to the parton saturation phenomena could 
be present [3] . Following this point of view, recently the structure function measured at HERA 
has been studied considering an analysis based on the first non-linear corrections to the DGLAP 
evolution [4], showing for the gluon distribution the nonlinear effects play an increasingly role 
at x < 10~ 3 and Q 2 < 10 GeV 2 . The procedure gives stable results for that distribution at both 
low x and Q 2 , instead of negative or constant values as in the recent NLO DGLAP analysis 
(see [4] for details). 

On the other hand, at very small x the leading logarithms a s ln(l/x) are shown to be 
important. In the leading logarithmic approximation (LLA) the QCD Pomeron corresponds 
to the sum of ladder diagrams with reggeized gluons along the chain, which are strongly or- 
dered in momentum fraction x. Such sum is described by the Balitzkij-Fadin-Kuraev-Lipatov 
(BFKL) equation [5]. The corrections at next-to-leading level (NLLA) are now known, lead- 
ing to strong modifications in the LLA spectrum [6]. At first glance, the results turned out 
conflicting and not useful for suitable phenomenological studies. Further, these features have 
been cured by canceling not desired singularities via a suitable resummation at all orders of 
the perturbative expansion, satisfying the renormalization group requirements [7, 8, 9]. The 
main shortcoming in doing phenomenology with a BFKL evolution beyond the leading order 
is the complicated mismatch between QCD perturbative and non-perturbative inputs, as the 
corresponding factorization properties are more involved than the DGLAP evolution. 

A very promising approach encoding all order a s ln(l/x) resummation is provided by the 
QCD dipole picture. It was proven that such approach reproduces the BFKL evolution [10]. 
The main process is the onium-onium scattering, that is the reaction between two heavy quark- 
antiquark states, in such way that the process is perturbative since the onium radius provides 
the essential scale at which the strong coupling constant is evaluated. In the large N c limit, 
the heavy pair and the soft gluons are represented as a collection of color dipoles. The cross 
section is written as a convolution between the number of the dipoles in each onium state 
and the basic cross section for dipole-dipole scattering due to two-gluon exchange. The QCD 
dipole model can be applied to deep inelastic lepton-nucleon scattering, assuming that the 
virtual photon at high Q 2 can be described by an onium. On the other hand, the proton is 
approximately described by a collection of onia with an average onium radius to be determined 
from phenomenology. This model has produced a successful description of the old structure 
function data [11]. 

The QCD dipole picture allows a systematic framework for testing the resummed next-to- 
leading order BFKL evolution kernels, producing predictions for the proton structure function. 
A method for doing this was recently proposed [12], where the resummation schemes are tested 
through the Mellin transformed j-moments of F2. Moreover, it has been shown that a geometric 
scaling for the photon-proton cross section and the symmetry between low and high Q 2 regions 
are associated to the symmetry of the two-gluon dipole-dipole cross section [13]. There, the 
proton is assumed to be a collection of independent dipoles at the time the interaction whose 
sizes are distributed around 1/Q s (x), with the latter being also the mean distance between the 
centers of the neighboring dipole. The quantity Q 2 (x) = A 2 e A log ( 1 / x ) is the saturation scale 
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[14]. Further, it was shown that the local geometric scale can be testable experimentally [15]. 
That is, the cross sections for exclusive processes only depend on the ratio of scales Q/Q s (x, b), 
where now the saturation scale Q 2 (x,b) has a impact parameter (6) dependence. 

In this work we revisit the phenomenology using this successful approach. The analysis is 
organized as follows. In the next section we shortly review the main formulae for the QCD 
dipole picture and give the phenomenological expressions for the structure functions and for the 
gluon distribution. In Sec. 3 we analyze the fitting procedure to the recent structure function 
data, determining the quality of the fit procedure and the range of applicability of the model. 
At the last section we drawn the main conclusions and summarize the results. 



2 The QCD color dipole picture and structure functions 



Let us start by describing the main features and basic formulas of the QCD dipole picture at 
high energies. The forward amplitude, A (s, t), for the onium-onium scaterring, integrated over 
impact parameter is written as, 



A (s, t = 0) 



/ ld 2 ri d 2 r 2 J dz 1 dz 2 $W(r 1 , Zl ) $ (0) (r 2 , z 2 ) 
d 2 s 1 d 2 s 2 



x 



Sl s 2 



n ( Y / 2 , r i, s i) n 07 2 , r 2 , s 2 ) a dd (s u s 2 ) , 



(1) 



where $^(rj, Zj) is the squared wave function of the quark-antiquark part of the onium wave 
function, r*j being the transverse size of the quark-antiquark pair and Zi the longitudinal mo- 
mentum fraction of the antiquark. The rapidity is given by Y/2 ~ ln(x /x), where x is a 
phenomenological constant and x is the Bjorken variable labelling the softer end of the pro- 
duced dipole. The dipole density is labeled as n {Y/2, r i: Sj), which is discussed below. The 
quantity Odd is the elementary dipole-dipole cross-section and reads as, 



Odd (*i,*2) = 2 7ca 2 [min (s!,s 2 )] 



In 



/ max (si, s 2 ) 
\ min (si, s 2 ) 



(2) 



In the large N c limit and in the leading-logarithmic-approximation the radiative corrections 
are generated by emission of gluons with strongly ordered longitudinal momenta fractions » 
Zi+i. The onium wave function with n soft gluons can be calculated using perturbative QCD. 
In the Coulomb gauge the soft radiation can be viewed as a cascade of color dipoles emitted 
from the initial dipole, since each gluon acts like a quark-antiquark pair. The dipole density 
n(Y, r, s) was defined [10] in such way, 

N (Y, s) = JdzJ d 2 r $ (0) (r, z) n(Y, r, s) (3) 

is the number of dipoles of transverse size s with the smallest light-cone momentum in the pair 
greater than or equal to e~ ¥ p + , where p + is the light-cone momentum of the onium. The whole 
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dipole cascade can be constructed from a repeated action of a kernel /C on the initial density 
n (r, s) through the dipole evolution equation, 



rY poo 

n(Y,r,8 1 )=n (r,8i)+ dy ds 2 JC(s 1: 8 2 )n(y,r, s 2 ) . (4) 

JO JO 

The evolution kernel /C is calculated in perturbative QCD. For fixed a s and in the limit 
of large N c the kernel has the same spectrum as the BFKL kernel [10]. Therefore, the two 
approaches lead to the same phenomenological results for inclusive observables. The solution 
of Eq. (4) is given by [10], 

„ (Y r s) - - r ex P(""' r > exD ( h*' (»•/») ^ l5) 

where ujp = 4 a s In 2, with ci s = a s N c / ir, and ajp = 1 + cjp is the Pomeron intercept. Therefore, 
the cross section grows rapidly with the energy since the number of dipoles in the light cone wave 
function grows rapidly on energy. The result in Eq. (5) has been obtained considering a process 
having only one scale, that is the onium radius. In the electron-proton deep inelastic scattering 
two scales are present and this result should be affected by non-perturbative contributions. The 
hard scale is given by the photon virtuality and the soft one is associated to the proton typical 
size. 

Now, let us review the obtention of the proton structure functions in the QCD dipole picture. 
In order to do so, an assumption is made that the proton can be approximately described by 
onium configurations. The basic assumption is given by, 

aZ P T (x, Q 2 ) = al — (x, Q 2 , k 2 ; /4) h (fc 2 ; /4) , (6) 

where n (k 2 ; n 2 F ) is the probability of finding an onium in the proton, with fj/p being a unknown 
perturbative scale characterizing the average onium size, as will be seen later on. We follow 
Refs. [11] to write down the virtual photon-onium cross section, where the /c^-factorization 
approach [16] was considered within the framework of the QCD dipole model. It reads as, 



,Q 2 ) = J d 2 r dz $ (0) (r, z) a^ dtpole (x, Q 2 , r) , (7) 



where (r, z) was already defined before. The virtual photon-dipole cross section is written 

as, 

Q 2 a'~ d ^ (x, Q 2 , r) = f d 2 k £ ^ a ig ( ^ T omum (z, k, r) . (8) 

where a ig /Q 2 is the Born cross section of an off-shell gluon of transverse momentum k in the 
the subprocess 7*g — > qq and T on i um is the unintegrated gluon distribution function of an onium 
state of size r, containing the hard Pomeron physics. The relation between this quantity and 
the dipole density is written as, 

k 2 T omum (z, k,r)= j ^ £ ^ n (z', r, a) a^ ole , {a fc) 2 ) , (9) 
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where n (z f , r, s) is the density of dipoles of transverse size s with the smallest light-cone 
momentum in the pair equal to z'p + in a dipole of transverse size r, of total momentum p + . 
This is given by the solution in Eq. (5). 

The next step is to apply a Mellin transform in the x-space and further a second Mellin 
transform in the A; 2 -space, where the conjugated variables are to and 7, respectively. After that, 
one obtains [11], 

Q 2 a^ ole (x,Q 2 ,r) = j p- j p- o ig (u, 1 - 7) Fonium (w, 7) (Q 2 r 2 ) 1 x~" . (10) 



The unintegrated gluon distribution for an onium of size r can be obtained from perturbative 
QCD and it is written as [10], 

T ( \-~ 1 ni x 

onium yd ; 7/ _ / \ ' \ ) 

7 w - a s x{v 

where in the leading-order BFKL approach, Xhoil) = 2^(1) — ip (7) — xjj (1 — 7) and a s x (7) 
is the eigenvalue of the LO BFKL kernel. The Mellin transform of the gluon-dipole coupling, 
v (7) , can be obtained from /c T -factorization approach as, 

2 -2 7 -l Wl-7) 



The unintegrated gluon distribution, Eq. (11), has a pole at up = a s x(l) an d then the 
integration over oj can be done analytically. Furthermore, the cross section be written in 

function of the LO photon impact factors /il,t(7), corresponding to the perturbative coupling 
to the photon. Namely, a ig = 4 n 2 a em e 2 h L T (7), with e 2 being the total charge of the quarks 
flavor contributing to the reaction and where, 



h T \ _ a s [r(i- 7 )r(i + 7 )] 3 1 /(i + 7 )(i-|) 



h L J 97 r (2-27) r (2 + 27) 1- 1 7 V 7(1-7) 



(13) 



Putting Eqs. (13) and (11) in Eq. (10) and solving analytically on uj for the pole on u = Up, 
one obtains, 



Q 2 a^° le (x, Q 2 , r) = 8n 2 a em e) a s [ p, h L , T (7) ^ (r 2 Q 2 Y exp 

J J 2, Til 7 



a s X (7) In — 

x 



.(14) 



The virtual photon-onium cross section is obtained replacing Eq. (14) in Eq. (7). The 
result depends on the squared wave function <3>(°) of the onium state, which cannot be computed 
perturbatively. This is solved by eliminating that dependence through an averaging over the 
wave function of transverse size [11], 

(/4)-t = ( r 2 ) = J dzd 2 r (r 2 ) 7 $ (0) (r, z) , (15) 
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where /j,p is a scale which is assumed to be perturbative. With this result, the photon-onium 
cross section, Eq. (7), takes the following form, 



7 onium 



(x,Q 2 ) 



An 2 a e 
0^ 



2a. 



I 



dl , f , v (7) 

TT~. h L,T 7 

2 m 7 



/4 



cxp 



a s X (7) In 



x 



(16) 



In order to compute the photon-proton cross section one needs to know the probability of 
finding an onium in it, n (7; fi F ). Relying on renormalization group properties, namely the 
overall result should be independent of the scale perturbative pip, an suitable ansatz is given 
by [11], 



/4 



n (7; up) = n cS (7) ^ ^ 



(17) 



where Qq is interpreted as a characteristic non-perturbative scale of the hadron-onium scat- 
tering. The quantity n e fr means the average number of primary dipoles in the proton and 
r = 2/Q their average transverse diameter (see second reference on [11]). Placing the ansatz 
Eq. (17) in Eq. (6) and using the result on Eq. (16), the proton structure functions can be 
calculated. The integration over 7 can be performed using the steepest descent method, with 
the saddle point given by x' (is) = — ln(Q 2 / 'Ql) / 'ot s hi(xo/x). Using the expansion of the BFKL 
kernel near 7 = |, one obtains, 



Is - 2 



1 



(18) 



where k(x) = [a s 7((3) In ^] 1 is the diffusion coefficient at rapidity Y = ln(xo/x). Putting 
all together, the structure functions are written in the simple form [11], 



F TtL (x,Q 2 ) = H TtL 



a s 7r 3 ej n eS 



96 



X °Y ,P Q / o t \l 



2 11 Q 



(19) 



where Ht = 9/2 and Hl = 1, defining n e s = n e s (j s ) and r = r (7 s ). In the nest section 
we consider the equation above to determine the parameters of the model by fitting the recent 
data on the proton structure function F2 — Ft + Fl. The conditions to obtain the structure 
functions from the saddle-point method fix the conditions of its safe applicability. Hence, the 
following bound should be covered, 



k (x) In 




Mi 



(20) 



which is realized for the region of moderate Q/Qo when compared to the range on xq/x. It is 
worth to mention that at higher Q 2 , the DLLA becomes valid, with the pole now at I/7 in the 
kernel x (7)- 
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PAR A METER 


T 

1 


TT 

11 


TTT 

111 




\f 
j\p 




ft D7R7 


U.Uc/OO 


OPT) rlinnlp 


n„ 

VO 




ft t\AA 


U.JO / 




Xn 


1.90 


1.0 (fixed) 


1.0 (fixed) 




1 + uJjp 


1.25 


1.26 


1.23 




Ar 


1.10 


1.23 






dR 


-0.198 


-1.01 




Non-singlet 


d 


-1.16 


-0.416 






b 


-0.482 


-0.007 




xVd.o.f. 




1.26 


1.25 


1.02 



Table 1: Parameters of the fit for the HI data analysis. 

3 Fitting results and discussions 

In this section we will perform a fitting procedure using the recent HERA data on the proton 
structure function [17, 18] using the QCD dipole picture given by Eq. (19). Having determined 
the parameters of the model from data, the longitudinal structure function and the proton 
gluon distribution are calculated without further adjustments. We made an simple attempt to 
describe phenomenologically the charm structure function, although the full expressions for it 
to be known in the literature (last reference in [11]). 

For sake of simplicity, we have defined the overall normalization for the proton structure 
function, Af p = (Ht + Hi) a s ii z e 2 n e fr/96. Also, for some fitting procedures we will include a 
non-singlet contribution and threshold (large x) effects. Namely, the hard Pomeron expression, 
Eq. (19), is multiplied by the factor (1 — x) Us , while the non- Pomeron contribution is multiplied 
by (1 — x) n " s . The expressions for the non-singlet piece and n s ^ ns read as, 

F™{x,Q 2 ) = Arx 1 "^ [ ^ ) (l- x ) n ™ iQ2) (21) 




"•«« 2 » - + n -^-\( l + WT-i' (22) 

where ctR = 0.4 is the Reggeon intercept, considered fixed in this analysis. Other values can be 
considered, however the x 2 quality of the fit is worst. 

Here, we have considered the following fitting procedures: 



• (I) Adjusting data in the whole range 1.5 < Q 2 < 150 GeV 2 and all x, considering the 
QCD dipole, Eq. (19), plus the non-singlet contribution, Eq. (21). The motivation is to 
include the large x data at high Q 2 , where the fixed target data (E665 and NMC points 
[19, 20]) were also added. The results for HI and ZEUS data sets are presented in Tables 
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PAR A METER 


T 

1 


TT 

11 


TTT 

111 




\f 
j\p 


U.UOO 1 


u.uuoo 


ft DQ77 


OPT) rlinnlp 


n„ 

VO 


P) zL7^ 


fl ^1 7 


yj.o i i 




Xn 


1.30 


1.0 (fixed) 


1.0 (fixed) 




1 + uJjp 


1.28 


1.29 


1.24 




Ar 


1.42 


1.43 






dR 


0.383 


-0.273 




Non-singlet 


d 


-0.842 


-0.808 






b 


1.36 


-0.677 




xVd.o.f. 




1.27 


1.25 


1.08 



Table 2: Parameters of the fit for the ZEUS data analysis. 

(1) and (2). The number of points is equal to 230 (HI) and 259 (ZEUS), respectively. 
The total number of parameters is 7; 

• (II) Adjusting data in the range 1.5 < Q 2 < 60 GeV 2 and all x, considering the QCD 
dipole plus the non-singlet contribution. We have fixed the parameter xq = 1 in this case, 
reducing the final number of parameters. The results for HI and ZEUS data sets are 
presented in Tables (1) and (2); 

• (III) Adjusting data in the whole range 1.5 < Q 2 < 150 GeV 2 and restricting to small 
region x < 10 -2 , considering only the QCD dipole contribution. We have considered the 
fixed values xq = 1 and n s = 7. 

Lets start discussing the procedure (I). The main difference between the two data sets are 
the values for the parameter x , being larger for ZEUS than for HI. On the other hand, the 
effective intercept for the hard QCD contributions is quite the same, that is 1 + oop = 1.25 (8), 
with the quality of fit better for ZEUS data set. It is worth to mention that the low values 
obtained for the parameters b and d suggests that they can be fixed in the usual way, n s = 7 
and n ns = 3, in the range considered here. This procedure would reduce strongly the number 
of free parameters. From the fit, it is also clear the non-singlet piece is Q 2 independent, since 
cir is quite small. 

For the procedure (II), the picture is similar to the previous one, with the effective intercept 
1 + ujjp = 1.26 (9) coming out different for the two data set. The quality of fit is also slightly 
equal in the two cases (see Tables). The conclusion about n s ns is the same as the procedure 

(I)- 

Finally, in the procedure (III) we have considered only the small x < 10 -2 data (102 data 
points for HI and 107 for ZEUS). In this case, xq = 1 and n s = 7 are taken as fixed, reducing 
the number of parameters equal to 3. Indeed, their values are in large extent similar with a 
better quality of fit for the HI data set. This procedure is similar to the previous analysis on 
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Refs. [11], with an even lower effective power (1 + ujjp ~ 1.282 in [11]). The low value for 
the fixed coupling constant a s ~ 0.1 reveals the well known necessity of sizeable higher order 
corrections to the approach. Accurate analysis in this lines, considering resummed NLO BFKL 
kernels, has been proposed recently [12] producing a reasonable a s ~ 0.2 for the typical Q 2 
range considered in phenomenology for structure functions. 

In order to illustrate the analysis, in Figs. 1 and 2 we present the fit result (procedure II) for 
the inclusive structure function for HI [17] and ZEUS [18] data set, including the fixed target 
points [19, 20] and extrapolated up to Q 2 = 150 GeV 2 . It is worth to mention that the new data 
are more precise than the previous analysis and the same quality of fit is obtained in a small 
range on Q 2 , in contrast with the studies in [11]. However, a rough data description can be 
obtained up to higher values. Here, our aim is discuss the model as the quality of fit. We have 
checked using the new HI data that a xVdof = 1.4 is obtained in the range 1.5 < Q 2 < 150 for 
all x. This in some extent shows the range on virtuality where a LO BFKL approach is suitable. 
That is, low x and not so high Q 2 . In virtualities of hundreds of GeV, the physics should be 
described by the DLA limit of the BFKL approach and which is a common limit also in the 
DGLAP evolution [21]. An interesting further study is to realize a best fit test to the DLA limit, 
determining its range of applicability. After this discussion, having determined the parameters 
of the model from data on F 2 , it is possible determine without further adjustments the structure 
function Fl and the gluon distribution x G(x,Q 2 ). From Eq. (19) and the definition of the 
overall normalization, we can see that Fl = (2/11) F 2 , since Ht + Hl — 11/2 and Hl = 1. In 
Fig. (3) this result is shown using parameters from procedure II from the HI data set against 
the recent F L data [17]. The results are in good agreement with the recent data, obtained 
without further adjustments. 

We made an attempt to describe phenomenologically the charm structure function, despite 
explicit expressions for it to be available. One suppose the following relation between the F 2 
structure function and F| c , 



where we have used m c = 1.5 GeV and the factor 2/5 corresponds to e 2 /{J2 6/ + e 2 ) . The 
result is shown in Fig. (4), where the parameters from procedure II (ZEUS data set) were 
taken into account. The simple ansatz is reasonable only at low Q 2 , overestimating the data on 
higher virtualities. It should be stressed that this function can be calculated using the complete 
expressions for the photon impact factors, Ht,l (7; m f), which consider also the dependence on 
the quark masses, to/ (see last reference in [11]). 

As a final analysis, in the QCD dipole approach the gluon distribution function can be 
calculated in a straightforward way from F 2 . The result is also independent of the overall 
normalization, which contains part of the non-perturbative inputs of the model. The gluon 
function is given by, 




(23) 



xG(x,Q 2 ) 



[h T ( 1 = ls ) + h L (l = lsT 1 F 2 (x,Q 2 ), 
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97s l-§7s r(2-2 7s )r(2 + 2 7s ) 

a, i + f 7s -f 7 2 [r(i- 7s )r(i + 7s )] 3 



F 2 (x,Q 2 ). (24) 



In Fig. (5) we shown the gluon distribution function obtained using Eq. (24) and the pa- 
rameters of the procedure (III) for the HI and ZEUS data sets. The curves represent xG(x, Q 2 ) 
as a function of x for distinct virtualities, 1.5, 15 and 150 GeV 2 . It would be timely to compare 
those curves with the ones from using a resummed kernel. In order to compare these predic- 
tions with NLO QCD DGLAP fit results, in Fig. (6) one plots the QCD dipole predictions with 
the recent QCD analysis. In the upper plot, we show the results for the HI QCD fit (central 
value, without bands) in Ref. [22] for Q 2 = 4 GeV 2 and the earlier analysis in Ref. [23]. The 
deviation between the recent HI fit and the QCD dipole curves (parameters from procedure 
III) is sizeable, whereas it is consistent with the previous analysis. Deviations are found in the 
comparison with the ZEUS NLO QCD fit [24] at Q 2 = 7 GeV 2 , in the lower plot. However, it 
should be stressed that the gluon distribution is an indirect observable, and its determination is 
model-dependent. For sake of comparison with other determinations of the gluon distribution, 
we have plotted the result using the Regge-like hard Pomeron contribution [25], which it is 
shown in the lower plot for Q 2 = 7 GeV 2 . Such approach shows that the ratio of the gluon 
distribution to the hard- Pomeron part of the singlet quark distribution is about 8 ± 1 in a large 
range on virtualities. Furthermore, the dependence on Q 2 of the hard-Pomeron contribution 
presents very good agreement with DGLAP evolution for Q 2 > 5 GeV 2 . 



4 Conclusions 



We have revisited the phenomenology to the structure functions in the QCD dipole picture, 
performing a fit to the recent data on the proton structure function. The model for the hard 
Pomeron provides a reduced number of free-parameters and the results produces a good quality 
fit. The adjusted constants are the overall normalization to F 2 , the effective Pomeron intercept 
ap = 1 + oj]p, the average transverse diameter of the primary dipoles in the proton r = 2/Qo 
and Xq scaling the ln(l/x) behavior. Moreover, the longitudinal structure function and the 
gluon distribution function can be determined without further adjustments. The results for Fl 
are in good agreement with the recent measurements. The main shortcoming from the approach 
is a low value for the strong coupling constant, a s ~ Oil, showing that a next-to-leading level 
analysis must be performed. The number of initial dipoles is either small, n e // = 3, and they 
have a large size r — 1 fm. Resummed perturbative kernels should provide a consistent value 
for the strong coupling constant in the range of virtualities considered, as reported in Ref. [12]. 
In this sense, the model presented here, although providing a remarkable data description, 
can be considered as effective and further phenomenology should consider higher order effects. 
Regarding the gluon distribution, the result is distinct from the recent NLO QCD fits, producing 
higher values at small x. We have also made a simple ansatz for describe the charm structure 
function, obtaining a very rough result. This is due we are disregarding the dependence of the 
photon impact factors on the quark mass, which modifies the final result in contrast with the 
assumption of avoiding the quark mass. 
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Figure 1: The inclusive structure function F 2 for the HI data set [17] and fixed target data 
points [19, 20]. The curve corresponds to procedure (II). 
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Figure 2: The inclusive structure function F 2 for the ZEUS data set [18] and fixed target data 
points [19, 20]. The curve corresponds to procedure (II). 
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Figure 3: The longitudinal structure function Fl [17], where the parameters of procedure II 
(HI data set) were taken into account and considering Fl = (2/11) F 2 . 
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Figure 4: The charm structure function F| c , where the parameters of procedure II (ZEUS data 
set) were taken into account. 
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Figure 5: The gluon distribution function x G(x,Q 2 ), where the parameters of procedure III 
were taken into account. In the plot on the left, parameters from HI data set fit and on the 
right the ZEUS data set fit. 



17 




Figure 6: The gluon distribution function xG(x,Q 2 ): comparison among HI and ZEUS QCD 
DGLAP fit and the QCD dipole prediction. 
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